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FOREWORD 


This  is  an  interim  report  on  the  study  of  underwater  noise  reduction 
techniques  'being  conducted  by  the  Sperry  Division  for  the  Office  of  Naval  Research, 
Acoustics  Programs,  under  Contract  No.  N00014-67-C-0303.  Included  in  this  report 
is  a  description  of  the  work  being  done  on  various  types  of  acoustic  baffles.  The 
theoretical  analysis  of  a  multilayer  cylindrical  baffle  is  given;  spring- mass  type  baffles 
are  analyzed  using  a  lumped-parameter  approach.  For  both  of  these  baffles,  good 
agreement  is  obtained  between  analytical  and  experimental  data.  Also  presented  are 
both  theoretical  and  experimental  results  of  a  new  baffle  material,  Sonite,  which  shows 
promise  of  operating  at  deep  depths  without  the  need  for  pressure  compensation. 
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SECTION  1 


INTRODUCTION 


The  effectiveness  of  the  acoustic  baffle  is  often  a  determining  factor  in  the 
maximum  capability  of  a  sonar  system.  Basically,  the  baffle  is  a  structure,  placed 
between  an  unwanted  signal  or  noise  source  and  a  sensor,  to  shield  the  sensor  from  the 
unwanted  source.  To  achieve  an  array  gain  consistent  with  the  directivity  index  of  the 
array,  efficient  baffle  structures  must  be  provided.  A  further  purpose  of  the  baffle  is 
to  provide  a  desired  response  characteristic  for  the  hydrophone  array.  Thus,  the 
acoustic  baffle  is  a  primary  tool  for  underwater  noise  reduction  and  directivity  pattern 
control  in  sonar  arrays. 

The  transmission  of  acoustic  energy  in  water  can  be  impeded  by  a  compliant 
baffle.  To  be  effective,  such  a  baffle  must  be  sufficiently  stiff  to  withstand  hydrostatic 
pressures  yet  it  should  be  dynamically  soft  to  perform  as  a  highly  reflective  structure. 

A  reflective  characteristic  is  acceptable  for  the  side  of  the  baffle  facing  the 
interfering  noise  source,  but  the  portion  facing  the  hydrophone  presents  a  major  design 
challenge.  Unless  this  side  is  nonreflective,  the  interaction  between  the  incident  signal 
and  the  reflected  signal  from  the  baffle  can  alter  the  hydrophone  array's  response,  often 
seriously  degrading  it. 

Under  the  current  program,  the  Sperry  Division  has  developed  a  set  of 
analytical  models  for  predicting  the  reflectivity,  transmittivity,  and  other  acoustic 
parameters  of  planar  multilayered  baffle  structures.  This  capability  extends  to  the 
prediction  of  directivity  pattern  and  response  for  line  hydrophones  at  any  steering  angle, 
including  end  fire. 

The  analytical  models  have  been  incorporated  into  a  FORTRAN  IV  computer 
program  which  has  been  described,  together  with  its  implementation  and  built-in  plotting 
routines,  in  previous  interim  reports  (references  1  and  2).  Empirical  and  theoretical 
data  have  been  successfully  correlated  with  the  results  of  the  computer  trials. 
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Several  aspects  of  baffle  technology  are  discussed  in  this  report.  First,  the 
analysis  of  a  multilayer  planar  baffle  was  extended  to  a  multilayer  cylindrical  baffle. 

As  in  the  planar  case,  the  cylindrical  baffle  consists  of  random  layers  of  fluid  and 
elastic  media.  Secondly,  baffle  designs  based  on  the  use  of  mechanical  springs  coupled 
to  masses,  which  represent  a  compliant  steel  plate  structure,  were  also  investigated. 
Lumped  parameter  models  were  developed,  which  were  used  to  calculate  insertion  loss. 
Results  of  studies  on  a  new  particulate,  fiber-reinforced,  ceramic  baffle  material  are 
presented  in  Section  4.  This  material  can  be  fabricated  into  intricate  shapes  with  con¬ 
ventional  techniques  and  shows  promise  of  operating  as  an  effective  baffle  at  high  pre- 
sures  without  requiring  pressure  compensation. 


SECTION  2 


MULTILAYER  CYLINDRICAL  BAFFLES 


2.  1  INTRODUCTION 

The  analysis  of  cylindrical  baffles  is  a  natural  extension  of  the  work  done 
on  planar  baffles  during  previous  phases  of  this  program.  In  this  section  an  infinitely 
long  cylinder,  with  wall  composed  of  layers  of  various  materials,  is  taken  to  represent 
the  baffle  configuration.  A  plane  wave,  parallel  to  the  cylindrical  axis,  is  the  excita¬ 
tion  (see  figure  2-1).  An  analysis  was  performed  to  determine  the  acoustic  pressure  at 
any  point  outside  of  the  cylinder.  With  this  information,  the  response  of  a  hydrophone 
placed  outside  the  baffle,  as  well  as  the  transmission  loss  through  the  baffle,  can  be 
computed. 


Solution  of  the  wave  equation,  consistent  with  the  boundary  conditions  im¬ 
posed  by  the  multilayer  cylinder,  forms  the  basis  of  the  analysis.  The  fundamental 
equations  for  this  approach  are  contained  in  Section  2.2  of  reference  3. 

A  description  of  a  computer  program,  which  calculates  the  response  of  a 
hydrophone  outside  the  cylindrical  baffle,  is  presented  in  this  section  following  the 
theoretical  analysis.  The  results  of  the  program  are  compared  to  experimental  mea¬ 
surements  and  arc  also  presented.  A  list  of  symbols  is  included  at  the  end  of  the  sec¬ 
tion  (paragraph  2. 7). 

2.2  BASIC  SOLUTIONS  TO  THE  WAVE  EQUATION 

The  equation  of  motion  for  a  homogeneous  isotropic  elastic  medium  is 

pi i  =  (\  +  2n)  grad  div  u  -  g  curl  curlu  (1) 

where  u  is  the  particle  displacement  vector,  \  and  nare  constants  characteristic  of  the 
medium,  and  p  is  the  density.  It  can  be  shown  that  the  most  general  solution  of  (1)  has 
the  form 


u  =  grad  P  +  curl  Q 


(2) 
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where  the  scalar  P  and  the  vector  Q  satisfy 


(3.1) 

(3.2) 


pP  =  (V  +  2 )jl)  V2P 

pQ  =  m-V2Q 

In  the  case  of  a  fluid,  m-  is  set  equal  to  0  and  \  is  replaced  by  the  bulk 

modulus  k. 

For  the  present  application,  equations  (3)  may  be  specialized  in  two  ways. 
First,  motion  at  a  single  frequency  will  be  assumed;  the  components  of  u  will  then  be 
complex  numbers  representing  the  phase  and  amplitude  of  the  displacement  components 
along  the  coordinate  axes.  The  actual  physical  displacement  vector  is  the  real  part  of 
u  exp  (-i  wt).  Second,  it  will  be  assumed  that  all  components  of  u  are  independent  of  z, 
and  that  the  z-component  vanishes  identically.  These  are  natural  assumptions  when 
considering  normal  incidence  on  an  infinite  cylinder.  One  may  then  write 

Q  =  S  'z'  (4) 

curl  Q  =  grad  SX1^ 

where  z  is  the  unit  vector  in  the  z  direction.  Equations  (3)  now  take  the  form 

V2P  +  n2P  =  0  (5. 1) 

V2S  +  m2S  =  0  (5.2) 

where  the  compressional  and  rotational  wave  numbers  n  and  m  are  given  by 

n2  (\  +  2|i)  =  w2p  =  m2(j.  (6) 

The  multilayer  cylindrical  baffle,  with  the  z  coordinate  suppressed,  becomes 
in  effect  a  set  of  concentric  rings.  Solutions  to  equations  (5),  valid  for  specific  values 
of  m  and  n,  can  only  apply  to  the  zone  between  two  circular  boundaries.  To  facilitate 
matching  the  solutions  on  opposite  sides  of  a  boundary,  so  that  continuity  of  stress  and 
displacement  can  be  maintained,  the  use  of  polar  coordinates  is  indicated.  Solutions  of 
the  form  P  =  BR,  where  B  is  a  function  of  bearing  b  only,  and  R  depends  only  on  range 
r,  will  be  especially  useful.  Transforming  the  operator  v2  into  polar  coordinates, 
equation  (5. 1)  gives 

B"/r2B  +  R’/rR  +  R”/R  =  -n2  (7) 

as  the  condition  for  B  and  R. 
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Since  B"/B  in  (7)  can  be  expressed  as  a  function  of  r  only,  while  in  fact  it 
is  a  function  of  b  only,  it  must  be  a  constant.  Hence  B  must  be  a  linear  combination  of 
sin  pb  and  cos  pi)  for  some  number  p.  Furthermore,  when  the  bearing  b  increases  by 
2t t,  the  values  of  B  must  repeat.  Thus,  p  must  be  an  integer.  The  condition  for  R 
is  then 

R”  +  R'/r  +  [n2  -  (p/r)2]  R  =  0  (8) 

which  requires  R  to  be  a  linear  combination  of  Jp(nr)  and  Hp(nr);  by  Jp  is  denoted  the 
Bessel  function  of  the  first  kind,  and  by  Hp  the  Hankel  function  Hp(l). 

Replacing  n  by  m  provides  solutions  S  to  equation  (5.2).  The  next  step  is  to 
form  more  complex  solutions,  for  both  P  and  S,  by  combining  the  basic  solutions.  For 
instance,  letting  Bp(b)  denote  a  linear  combination  of  sin  pb  and  cos  pb,  and  using 
Rp(x)  to  indicate  a  linear  combination  of  Jp(x)  and  Hp(x),  the  series 

£  Bp(b)  Rp(kr)  (9) 

P=o 


will  provide  a  general  solution  for  P  when  k  =  n  and  for  S  when  k  =  m. 

A  solution  of  the  form  (9)  is  required  for  each  zone  of  constant  n  and  m. 

The  constants  used  in  forming  the  linear  combinations  Bp  and  Rp  must  be  adjusted  to 
make  the  solutions  match  at  the  boundaries  between  zones.  Section  2.  3  discusses  how 
the  matching  is  done. 

2.3  BOUNDARY  CONDITIONS 

At  the  boundaries  separating  regions  of  different  wave  velocities,  two  types 
of  conditions  must  be  met.  First,  the  stress  vector  must  be  continuous  across  the 
boundary.  Second,  the  displacement  normal  to  the  boundary  must  be  continuous  as  the 
boundary  is  crossed;  and  in  the  case  of  a  solid-solid  boundary,  the  tangential  displace¬ 
ment  must  also  be  continuous  (assuming  secure  bonding  of  the  surfaces). 

The  stress  vector  t  acting  on  a  surface  may  be  expressed  in  terms  of  the 
displacements  by  the  formula 

t  =  a  \  div  a  +  (i  a  •  (7u  +  u  V)  (10) 

where  a  Is  the  unit  normal  to  the  surface.  Using  polar  coordinates  and  setting  a  equal  to 
r  (the  unit  vector  in  the  range  direction),  the  stress  components  normal  and  tangential 
to  the  boundary  are 

t  =  v  div  u  +  2n  du  dr 
r  -  r 

tp  =  g  j^(l  r)  dur  d  b  -  ub  r  +  dub  dr 
where  ur  and  ub  are  the  displacements  normal  and  tangential  to  the  boundary,  respectively. 
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At  the  outer  boundary,  the  displacements  and  stresses  must  be  compatible 
with  the  assumed  input  conditions.  In  this  analysis,  an  incident  plane  wave  parallel  to 
the  cylinder  axis  is  assumed.  The  displacements  in  the  fluid  surrounding  the  cylinder 
are  thus  the  sum  of  an  incident  plane  wave  and  waves  scattered  from  the  cylinder. 
Assuming  that  axis  b=o  points  toward  the  source  of  the  incident  wave,  it  may  be  de¬ 
scribed  as 


Pqj  =  exp  (-inr  cos  b) 

=  £  «p  (-i)P  Jp  (nr)  cos  pb  (12) 

o 

where  <p  equals  2  for  positive  p  and  1  for  p=o.  The  scattered  waves  may  be  written 


00 

£  cp  Hp  (nr)  cos  Pb 
o 


(13) 


where  the  coefficients  Cp  are  to  be  determined.  Note  that  the  scattered  waves  do  not 
make  use  of  a  fully  general  expansion  in  the  form  (9).  The  function  Rp  is  restricted  to 
a  pure  Hankel  function  so  that,  for  large  r,  each  term  of  (13)  will  correspond  to  an 
outgoing  circular  wave.  The  function  Bp  is  restricted  to  the  cosine  for  symmetry 
reasons:  with  the  incident  wave  coming  in  along  the  b=o  axis,  displacement  and 
stresses  normal  to  the  circular  boundaries  will  be  even  functions  of  b,  and  displace¬ 
ments  and  stresses  tangent  to  the  boundaries  will  be  odd  functions  of  b.  The  relations 


ur  =  dP/dr  +  (1/r)  dS/db 
ub  =  (i/r)  dp/ab  -  as/ar 


(14) 


and  also  equation  (11),  indicate  that  the  functions  P  must  have  even  symmetry  in  b, 
while  the  functions  S  must  have  odd  symmetry.  In  short,  the  expansions  of  form  (9) 
for  P  and  S  have  predetermined  choices  of  Bp;  only  Rp  is  open  to  selection.  Except 
for  the  restriction  on  the  scattered  waves,  one  may  generally  choose  Rp  freely;  how¬ 
ever,  in  the  innermost  zone  one  must  select  pure  Bessel  functions  to  avoid  a  pole  at 
the  origin. 


The  entire  process  of  matching  displacements  and  stresses  at  the  boundaries 
must  be  carried  out  separately  for  each  p,  until  the  terms  in  (12)  become  negligible. 

For  a  given  p,  let  F  and  G  denote  the  choices  for  Rp  in  the  expansions  of  P  and  S, 
respectively;  then 


F(r)  =  Cn  Jp  (nr)  +  C12  Hp  (nr) 
G(r)  =  C21  Jp  (mr)  +  C22  Hp  (mr) 


(15) 


iPjP 

I 
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Substituting  P  =  F  cos  pb  and  S  =  G  sin  pb  into  (14)  and  (11),  and  noting  that  div  u  =  -n^p, 
gives  the  pth  terms  in  the  sine  or  cosine  series  for  the  displacement  and  stress~com- 
ponents  Using  capital  letters  to  denote  the  series  coefficients,  so  that,  for  instance 


ur  =  Y,  Urp  cos  pb 
o 

one  finds  directly  that 

Urp  =  F'  +  (P/r)  G 

ubp  =  -  [G'  +  (PA)  f] 


Some  additional  calculation  shows  that  the  coefficients  in  the  expansions  for 
the  stress  components  are 


Trp  =  ^  +  2'iF”  +  2a  (p/r2)  (rG’  -  G)  (17.1) 

rbp  r  -t*n)2G  -  2  G”  -  2  .  (p  r2)  (rF*  -  F)  (17.2) 

where  (8)  has  been  used  to  modify  equation  (17.2). 


It  is  the  quantities  given  in  (16)  and  (17)  which  must  show  continuity  at  the 
boundaries,  with  the  exception  of  U^p  in  the  fluid/solid  case.  The  next  step  is  to  repre 
sent  these  displacement  and  stress  coefficients  explicity  as  linear  functions  of  the  con¬ 
stants  Cy  appearing  in  equations  (15).  This  will  be  done  in  Section  2.  4,  where  the 
procedure  for  solving  the  equations  is  discussed. 

2.4  NUMERICAL  METHODS 


UP011  substituting  (15)  into  (16)  and  (17)  the  following  explict  formulas  are 
obtained  giving  the  Fourier  coefficients  of  the  displacement  and  stress  components. 

The  symbols  x  and  y  are  used  for  nr  and  mr,  respectively. 

rUrp  -  x  Jp'(x)  C„  +  x  Hp’(x)  Cj2  +  P  Jp(y)  C21  +  p  Hp(y)  C22  (18.1) 

-rUbp  =  p.Ip(x)  Cy  +  p  Hp(x)  C12  +  y  Jp’(y)  C2i  +  y  Hp’(y)  C22  (18.2) 

,  2  Frp  =  [2  .1p"(x)  -  Jp(x)  ]  C,  j 

+  x2  [2>‘Hp"(x)  -  VHp(x)]c12 

+  2  *p  [y  •ip’(y)  -  Jp(y)j  c2i 

*  2  P  [  y  Hp’(y)  -  Hp(y)  J  c22  (18.3) 
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-r2  Tbp  =  2  |aP  [  x  Jp'(x)  “  Jp(x)  ]  cn 
+  2  up  [x  Hp'(x)  -  Hp(x)  j  C12 
+  i^y2  [2  Jp"(y)  +  Jp(y)  J  C2J 

+  py2  [  2  Hp"(y)  +  Hp(y)]  C22  (18.4) 

As  pi’eviously  stated,  all  the  quantities  (18)  except  (18.2)  must  show  con¬ 
tinuity  at  every  boundary,  and  at  a  solid/solid  boundary  (18.2)  also  must  show  continuity. 
The  meaning  of  this  is  that  if  a  quantity  such  as  (18.  3)  is  evaluated  on  the  outer  side  of 
a  boundary  r  =  r0  using  the  physical  constants  \J,  pi,  n1,  ml,  x1  =  nlr0  and  yl  =  mlr0 
pertaining  to  the  outer  medium,  and  the  values  Cjjl  describing  the  wave  in  the  outer 
medium,  then  the  resultant  value  of  (18.3)  will  be  the  same  as  that  found  on  the  inner 
side  of  the  boundary  using  the  physical  constants  \H,  pH,  n11,  m11,  xll  =  nllr0  and 
yll  =  mllr0  pertaining  to  the  inner  medium  and  the  values  Cij*1  describing  the  wave  in 
the  inner  medium. 

Thus,  at  each  boundary  a  set  of  equations  may  be  set  up,  expressing  the 
requirements  of  continuity.  The  equations  may  be  written 

(r  Urp)!  -  (r  Urp)n  =  0 

ubp^  ”  ("r  Hbp)**  = 

(r2  Trp)1  -  (r2  Trp) 11  = 

(-^Tbp)1-^  r2Tbp)n 

where  superscripts  I  and  II  indicate  evaluation  in  the  outer  and  the  inner  medium, 
respectively.  It  must  be  noted,  however,  that  not  all  the  equations  are  required  at 
every  boundary.  At  a  fluid/fluid  boundary,  only  the  first  and  third  equations  are 
needed,  since  the  fourth  is  satisfied  automatically  when  p  =  o  in  both  media.  At  a 
fluid/solid  boundary  all  but  the  second  are  needed,  and  all  are  needed  at  a  solid/solid 
boundary.  The  number  of  equations  (19)  needed  to  express  the  continuity  requirements 
can  thus  be  two,  three,  or  four,  depending  on  the  type  of  boundary. 

When  the  continuity  equations  have  been  written  for  all  the  boundaries,  one 
has  a  number  of  simultaneous  linear  equations.  The  unknowns  are  the  various  Cjj 
which  describe  the  waves  within  the  various  zones  of  constant  velocity.  There  will 
always  be  just  enough  equations  to  allow  a  solution  to  be  found,  based  on  the  given 
values  of  Cji  in  the  outermost  zone.  When  the  equations  have  been  solved  for  the 
value  of  C12  in  the  outermost  zone,  all  this  being  for  a  particular  value  of  p,  the  pth 
coefficient  in  the  expansion  (13)  has  been  calculated.  When  sufficient  coefficients  have 
been  determined,  the  overall  displacement  potential  in  the  outermost  medium  may  be 
calculated  from 


0 

0 

=  0 


(19) 
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P  =  poi  +  P 


02 


whore  the  notations  P01  and  PQ2  refer  to  (12)  and  (13),  respectively. 

The  actual  rms  pressure  in  the  outermost  medium  equals  times 

0.707  |  P|. 

2.5  COMPUTER  PROGRAM 


This  subsection  contains  an  annotated  listing  of  CYLBAF,  a  program  to  cal¬ 
culate  the  response  of  a  multilayer  cylindrical  baffle  to  normally  incident  waves.  The 
program  language  employed  is  time  sharing  FORTRAN. 

The  principal  function  of  the  program  is  to  find  the  coefficients  in  the 
scattet  ed-wave  expansion,  equation  (13).  These  coefficients  form  a  one-dimensional 
array  C02  in  the  program:  one  value  for  each  azimuth  harmonic  used  in  the  incident- 
wave  expansion,  equation  (12).  (The  order  of  the  harmonic  is  indicated  by  N  in  the 
program,  instead  of  the  symbol  p  used  in  the  equations.) 

To  find  the  scattered-wave  coefficient  for  any  N,  it  is  necessary  to  set  up 
the  matrix  A  showing  the  combined  boundary  conditions  (18)  for  all  the  boundaries  pre¬ 
sent.  The  structure  of  the  matrix  for  a  baffle  of  three  layers  is  illustrated  in  figure 
2-2.  Each  boundary  is  assigned  a  number  (NCON  in  the  program)  of  rows,  involving  the 
parameters  of  two  adjacent  zones.  Blocks  marked  ”1”  involve  the  parameters  of  the 
zone  on  the  inner  side  of  the  boundary,  while  blocks  marked  "11"  involve  the  zone  on  the 
outer  side.  Each  zone  is  assigned  NCOE  columns.  Unmarked  blocks  are  zero-filled. 

I' oi  each  marked  block,  sufficient  Bessel  functions  are  computed  to  allow 
the  coefficients  of  equation  (18)  to  be  calculated.  Then,  depending  on  the  type  of 
boundary,  the  proper  coefficient  values  are  put  into  the  matrix.  When  the  entire  matrix 
has  been  filled  in  this  way,  it  is  solved  for  the  scattered-wave  coefficient. 

F  inally,  the  coefficients  are  stored  on  a  disc-file.  Also,  the  pressure  at 
various  points  of  interest  is  calculated  and  stored  on  a  separate  file.  Thus,  immediate 
results  are  available,  but  also,  the  stored  coefficients  may  be  used  to  obtain  further 
results  at  a  later  time. 
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ANNOTATED  LISTING  OF  COMPUTER  PROGRAM 


♦CYL3AF:  PROGRAM  T0  CALCULATE  RESPONSE  0F  MULTILAYER  CYLINDRICAL 
BAFFLE  TO  NORMALLY- 1 NC I  DENT  PLANE  WAVE.  CAN  HANDLE  UP  TO  5  LAYERS. 

DIMENSION  IDENT(4>*THC7>*D<7>*C(7>*B<7)*LA<7)*MIJ<7)*TMU(7)* 

4  NC0E(7)*NC0NC6)»XS9<6*2>*X(6»2)»YSQ<6*2>#Y<6*2)*flF(5*12*4>* 

4  BV(2*3*4)*8VC(3*4)»A<23*24)*C02(50> 

REAL  LA  *  MU 

COMPLEX  3VC*A*C02*C0 1 
EQUIVALENCE  (BV*3VC> 

*USE  0F  BV  FACILITATES  HANDLING  REAL  AND  IMAGINARY  PARTS  OF  BVC. 

CALL  0PENF(1*"CYLDAT") 

♦*'C  YLDAT"  CONTAINS  PHYSICAL  PARAMETERS  DESCRIBING  BAFFLE  AND  INCIDENT 
WAVE*  FIRST  READ  RUN  IDENTIFICATION  AND  NUMBER  OF  CONSTANT-VELOC I TY 
ZONES*  INCLUDING  CORE  AND  EXTERNAL  ZONES. 

READ( 1  *  500 )  IDENTjNZON 
5OOF0RMAT(4A3*2XI  1  ) 

N3DD=  NZ0N-1 

♦NBDD  IS  THE  NUMBER  OF  INTERZONE  BOUNDARIES.  NEXT  READ  PARAMETERS 
FOR  EACH  ZONE:  TH=  LAYER  THICKNESS  (DUMMY  VALUE  REQUIRED  FOR  CORE 
AND  EXTERNAL  ZONES)*  D=  DENSITY*  C=  COMPRESS  1 3NAL-WAVE  SPEED* 

B=  ROTATIONAL-WAVE  SPEED. 

DO  5  K=  1  * NZON 
READ(1*)  Tl *T2*T3*T4 
♦FREE-FIELD  FORMAT 
TH(K)=  Tl 
D(K)=  T2 
C(K)=  T3 
5B(K)=  T4 

♦READ  OUTER  RADIUS  OF  CYLINDER.  READ  INCIDENT-WAVE  FREQUENCY  AND 
CONVERT  TO  RADIANS/SEC. 

READC!*)  R,T1 
0MEGA=  FPI C2*T1 ) 

♦FOR  EACH  ZONE  COMPUTE  LAME  PARAMETERS  LA  AND  MU*  AND  NC3E=  NUMBER 
OF  COEFFICIENTS  NEEDED.  FOR  EACH  BOUNDARY  COMPUTE  SPECIAL  SYMBOLS 
X  AND  Y*  AND  NC0N=  NUMBER  OF  CONTINUITY  CONDITIONS  ( EQS  lg>. 

T2=  R 

NC0N<  NBDD  )  =  1 
DO  10  L=  1  *  NZON 
K=  NZON-L+1 

LACK) =  D(K)*CCK>**2  -  TMU(K)=  2+MU(K>=  D(K>*3CK)*+2 
IF(L.E9  •  1 )GOTO  10 
Tl  =  0KEGA  +  T2 

XSQ(K*2)=  <  X ( K  *  2 ) -  Tl /CCK+l ) )**2 
YS9  ( K*  2  )  =  <  Y  ( K*  2  )  =  Tl  /B  <!<•*•  1  >)  **2 
♦SECOND  SUBSCRIPT  =  2  FOR  OUTER  MEDIUM*  =  1  FOR  INNE. 
vrO(K,l)=  (  X(  K* 1  ) =  T1/C(X))**2 
1  >=  ( Y ( K  *  1  >  =  Tl /BCK) ) **2 
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NC0EC  K )  =  4 

IF(B<K).EQ.O)  NC0E  ( K )  =  2 
NC0N(K)=  NCON(K)  «•  NC0E(K)/2 

♦THIS  IS  INNER-MEDIUM  CONTRIBUTION  T0  NC0N •  NOW  MAKE  OUTER-MEDIUM 
CONTRIBUTION  FOR  NEXT  BOUNDARY*  AND  UPDATE  RADIUS. 

IFCK.EQ.l )  GOTO  10 
NC0N<K-1)=  NC0E(K)/2 
T2  =  T2-TH(K) 

IOC0NTINUE 

♦ADJUST  NC0E ( 1 )  FOR  LACK  OF  HANKEL  FUNCTIONS  IN  CORE  ZONE. 

NCOEC 1 )  =  NC3EC 1 >/2 

♦FOR  EACH  ORDER  N.  WILL  NOW  WRITE  MATRIX  A  FOR  CONTINUITY  CONDITIONS. 
FIRST  MAKE  SURE  UNSPECIFIED  A-VALUES  ARE  ZERO. 

DO  12  1=1*23 
DO  12  J  =  1.2 A 
12A( I .  J)  =  0 
COl=  i 

♦C01  IS  COEFFICIENT  IN  INCIDENT-WAVE  EXPANSION  ( EQ  12>. 

N=  -1 
1 SN=  N*1 
I9M=  -1 

10=  J0=  1 

♦10  AND  JO  INDEX  UPPER-LEFT  ELEMENT  OF  MATRIX  BLOCK.  II  WILL  INDEX 
BOUNDARIES.  JJ  ZONES. 

DO  70  1 1  =  l . N9  DD 
DO  65  LL= I .2 

♦LL=  1  FOR  INNER  MEDIUM.  =  2  FOR  OUTER  MEDIUM. 

JJ=  1 1 ♦LL- 1 
KK=  2*11 +LL-2 

♦KK  IS  INDEX  FOR  BESS EL-FUNCTION  VALUES. 

IF C N. GT »0 ) GOTO  30 

♦FOR  N=0  COMPUTE  3ESSEL  FUNCTIONS  OF  ORDER  -2  THROUGH  2. 

DO  20  K=3 . 5 

CALL  BESSCBFCK.KK. 1 ).3F(K.KK.2>.K-3.XC I I.LL) ) 

20CALL  BESS (BF C K. KK. 3 ).3F<K.KK»  4) »K-3.Y< I I.LL)) 

DO  25  K  =  l  .4 

BFC l . KK. K ) a  BFC5.KK.K) 

2SBF ( 2. KK. K ) =  -9FC4.KK.K) 

GOTO  40 

♦IF  N>0 .  SHIFT  OVER  OLD  BESSEL  FUNCTION  VALUES  AND  COMPUTE  NEW 
FUNCTIONS.  OF  ORDER  N*2. 

30 DO  35  K  =  l .4 
DO  35  L=1 .4 

353F ( K»  KK.L ) =  8F<K*1.KK.L> 

CALL  BESS <BF( 5 . KK. I  ? . SF<  5. KK. 2 ) »N+2.  X ( I I.LL)) 

CALL  BESS<BF(5.KK,3).BF<5.KK,4).N*2.Y< II.LL) > 
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40SIGN=  3-2*LL 

♦THIS  IS  THE  CENTRAL  (MINUS)  SIGN  IN  ESS  19.  NOW  F0RM  SIGNED  J-  AND 
H-FUNCT I 0NS #  TOGETHER  WITH  FIRST  AND  SECOND  DERIVITIVES#  FROM 
THE  J-  AND  Y-  FUNCTION  VALUES  BF. 

DO  50  K=2#  4»2 
DO  45  L= 1 .2 

♦MANKEL  FUNCTIONS  FIRST. 

1=  K+L-2 

BV(L#1.K)=  SIGN*3F(3»KK»I) 

8VCL#  2#K)  =  .  5  +  SIGN*CBFC2#KK#I )-8F(4#KK# I) ) 

45BVCL#  3. K )  =  .25+SIGN+(8F(i#KK#  I  )-2*8F (3#KK#  I )  +3FC  5#  KK#  I ) ) 

♦THEN  J-FUNCTIONS. 

DO  50  J=1 # 3 

50BVCC J#K-1 )=  BVC  1  #  J#K) 

♦NOV/  FILL  IN  A  MATRIX  BLOCK.  IBM  =  1  FOR  11=  LL=  1#  WHEN  H-FUNCTIONS 
FORBIDDEN.  AFTER  THAT  IBM  STAYS  =  2. 

IBM=  (IBM+3J/2 
DO  55  13= 1  *  IBM 

*18=  1  FOR  J-FUNCTI0N  COEFFICIENT#  =  2  FOR  H -FUNCTION  COEFFICIENT • 
J=  JO ♦ IB- 1 

♦COMPRESS  I ONAL- WAVE  COEFFICIENTS#  EQS  18.1  AND  18.3s 
A(IO#J)=  X< I I#LL)*3VC(2#I8) 

ACIO+1 # J)=  XSQC  I  I#LL)*(TMU(JJ)*BVCC3# IB ) -LA( JJ ) *B VC < 1 # IB ) ) 
IFCNC0NC1I ) .EQ. 2) GOTO 55 
♦COMPRESS I ONAL-V/AVE  COEFFICIENTS#  EQ  18.4! 

ACIO  +2  #  J )  =  TMUC JJ)*N*(X(II»LL)*BVC<2# IB ) -B VC C 1 # 13 >  > 

IF(NC0N( II) .EQ.3)G0T0  55 
♦COMPRESS  I ONAL- WA  VE  COEFFICIENTS#  EQ  18.2; 

A< 10  +  3 #  J ) =  N*B VC  < 1 # IB ) 

55C0NTINUE 

IF(B( JJ) .EQ .0 )G0T0  65 
♦NOW  ROTATIONAL-WAVE  COEFFICIENTS#  IF  ANY. 

DO  60  13=1# IBM 
J=  JO+IB-1 +IBM 
♦EQS  18.1  AND  18.3! 

A<IO#J)=  N*BVC ( 1 # 1 3+2 ) 

ACIO+1 #J)=  TMUC  JJ)*N*(Y(II#LL)*3VC(2#  IB+2)-BVCC 1  # IB  +  2) ) 
IF(NCONdl)  .EQ.2)G0T0  60 
♦EQ  18.4! 

A CIO +2# J) =  YSQC I I#LL)*CTMUCJJ)*BVCC3# I B +2 )  +  MUC JJ ) *3 VC ( 1 » IB+2  )  ) 
IFCNC0NCII) .E3.3)GOT0  60 
♦EQ  18.2t 

ACIO  +  3#  J)  =  YC 1 1 #LL)*3VC(2# IB  +  2) 

60C0NTINUE 

♦UPDATE  JO  IN  PASSING  FROM  INNER  T0  OUTER  MEDIUM. 

65J0=  JO  +  C2-LL)*NC0EC II ) 

7010=  10  +  NCONCII) 
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♦ALL  MATRIX  ENTRIES  HAVE  BEEN  COMPUTED.  INTERCHANGE  LAST  TWO 
COLUMNS  OF  A*  AND  MULTIPLY  (NEW)  LAST  COLUMN  BY  -COW 
KK=  IO-NCONCNaDD) 

LL=  IO-l 
DO  75  I=KK,LL 
91  =  A ( I , JO  -  1 ) 

A(l.JO-l>=  -CO  1 *A ( I #  JO -2 ) 

7  5A ( I , JO -2  )  =  9  1 

♦NOW  SOLVE  MATRIX  FOR  SCATTERED-WAVE  COEFFICIENT  C02  CEQ  13). 

CALL  S0LV(A,  IG-1 ,C02(N*1 ) ) 

♦UPDATE  COl  AND  TEST  FOR  EXHAUSTION  OF  INCIDENT-WAVE  EXPANSION. 

IF ( N • GT .0 ) GOTO  gO 
COW  2 

80C0 I  a  COWCMPLXCO.,-1  .  > 

IF(ABS(3F(3.2*N3DD,1 ) > +ABS( BF< 4* 2^NBDD, I ) ) ♦A3SC3F ( 5,2+NBDD, 1 )) 
i.GT* 1 E-5 )GOT0  15 

♦CALCULATION  0F  SCATTERED-WAVE  COEFFICIENTS  COMPLETE*  STORE  THEM 
FOR  FUTURE  USE. 

CALL  0PENF(2,"CYLC0E"> 

1=  NM 

WR I TE (2 »  50  5 )  I  DENT. I 
50 SF0RMA  T ( 4A3 ,  1 4,  "  COEFFICIENTS") 

DO  85  J=1 , I 
§  1 =  C02CJ) 

85WRI TE (2*510 )  91 
5 1 0 FORMAT (2E 1 7 • 1 0 ) 

♦IF  IMMEDIATE  PRESSURE  READINGS  ARE  DESIRED#  TAKE  LOCATIONS  FOR 
READINGS  FROM  "CYLLUC"  AND  PUT  VALUES  INTO  "CYLVAL". 

CALL  0PENF(3,"CYLL3C") 

CALL  OPENF ( 4# "CYLVAL"  ) 

FACTOR3  SQRT( . 5 ) *D( NZ0N) ♦  OMEGAS 
WRITE(4,515)  IDENT 
51  5F0RMAT ("VALUES  FOR  ",4A3) 

90READC3,,END=999)  GAP, ANGLE 
Tl*  ( R+GAP ) ♦OMEGA /C ( NZON) 

T2=  FPI (ANGLE/180) 

T3=  -T1*C0S(T2) 

♦INCIDENT-WAVE  POTENTIAL: 

91=  CMPLX(C0S(T3),SIN(T3)) 

♦ADD  TF.RMS  OF  SCA  TTERED- WAVE  POTENTIAL  TO  GET  TOTAL  POTENTIAL* 

DO  95  J“0  » N 

CALL  BESS (T3,T4,J,T1  ) 

9591=  91  +  CC2  (  J+- 1  )  +CMPLX  (  T3»T4)*C0S(J*T2) 

♦COMPUTE  AND  STORE  PRESSURE  VALUE. 

T1  =  FACTOR+CABSCv l  ) 

WR  I  TE ( 4, 520 )  GAP, ANGLE, T1 
520F0RMAT(F7 .2,F  6.0, El  4.6) 

♦GO  BACK  AND  READ  NEXT  LOCATION. 

GOTO  90 

♦DESCRIPTION  OF  SUBROUTINES: 

SUBROUTINE  3ESSCB 1 ,B2,K, X)  COMPUTES  BESSEL  FUNCTION  B1  OF  FIRST  KIND 
AND  32  OF  SECOND  KIND  FOR  ARGUMENT  X  AND  ORDER  K. 

SUBROUTINE  SOLV(A,K,X)  TAKES  THE  FIRST  K  ROWS  AND  K+l  COLUMNS  OF  A 
AND  APPLIES  ROW  OPERATIONS  TO  MAKE  A(I,J)=0  FOR  I>J  AND  A(I,I)*|. 
THEN  X=  A ( K, K+l ) . 

999END 
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2.6  CORRELATION  WITH  EXPERIMENTAL  DATA 


Figure  2-3  shows  a  comparison  between  the  computer  program  prediction 
and  experimental  results  on  a  two-layer  cylindrical  baffle  consisting  of  5/8-inch  thick 
steel  and  I/2-inch  thick  RAL.  Good  agreement  exists  between  2  and  8  kHz;  below 
2  kHz,  the  theory  deviates  from  the  experimental  results  but  this  is  believed  due  to 
limitations  in  certain  subroutines  of  the  program  which  will  be  improved  in  the  future. 

Figures  2-4  and  2-5  show  the  predicted  and  experimental  azimuthal 
responses  of  a  hydrophone  in  front  of  the  same  baffle  as  in  figure  2-3  for  frequencies 
of  2500  Hz  and  4500  Hz,  respectively.  The  correlation  is  excellent.  Additional  com¬ 
parisons  of  this  type  are  necessary,  however,  before  a  proper  evaluation  of  the  program 
can  be  made. 

2.7  SYMBOLS 


a 

b 


B,  Bp 
Cp 

Cn,  C12,  C21,  C22 


F 

G 

Hp  =  Hp(l) 

Jp 

k 

S  H 

rn 

n 

to 


P 

Pol  ?o2 

Q 

r 

£ 

R,  Rp 
P 


Unit  vector  normal  to  a  surface 
Bearing  coordinate  in  polar  system  r,  b,  z 
Function  of  b  only 

Coefficient  in  scattered-wave  expansion 

Coefficients  in  expansion  of  F  and  G  (equation  15) 

Choice  of  Rp  for  compressional  waves 

Choice  of  Rp  for  rotational  waves 

Hankel  function  of  order  p 

Bessel  function  of  first  kind  and  order  p 

Bulk  modulus 

Lame  elastic  constants 

Rotational  wave  number 

Compressional  wave  number 

Angular  frequency  of  sound  wave 

Compressional-wave  potential 

Incident  and  scattered  wave  potentials 

Vector  potential  for  rotational  waves 

Range  coordinate  in  polar  system  r,  b,  z 

Unit  vector  in  range  direction 

Function  of  r  only 

Density 
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tSJ  <N| 


s 

t 

t 

V  lb 

Trp>  Tbp 
u 

“r.  ub 

^rp»  ^bp 
x  =  nr 

y  =  mr 


Scalar  potential  for  rotational  waves 
Time 

Stress  vector 

Stress  components  in  radial  and  tangential  directions 
Fourier  coefficients  of  tj.  and  t^ 

Displacement  vector 

Displacement  components  in  radial  and  tangential  directions 
Fourier  coefficients  of  ur  and 

Axial  coordinate  in  polar  system  r,  b,  z 
Unit  vector  in  axial  direction 
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Figure  2-1.  Structure  of  Three- Layer  Cylindrical  Baffle 


Figure  2-3.  Comparison  of  Predicted  and  Experimental  Results 
of  Two- Layer  Cylindrical  Baffle 
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FREQUENCY— KHz 


f  =  2500  Hi 

baffle 


THEORY 

EXPERIMENT 


STEEL  OUTSIDE 
RAL  INSIDE 

"  HYOROPHONE  SPACING 


Figure  2-4.  Predicted  and  Experimental  Hydrophone/Baffle 
Azimuthal  Response  at  2500  Hz 


f =4500  Hz 

baffle 


RADIUS  =  48" 

%"  STEEL  OUTSIDE 
V2"  RAL  INSIDE 
I  '/g”  HYDROPHONE  SPACING 


THEORY 

EXPERIMENT 


Figure  2-5.  Predicted  and  Experimental  Hydrophone/Baffle 
Azimuthal  Response  at  4500  Hz 


SECTION  3 


SPRING-MASS  BAFFLES 


3.  1  INTRODUCTION 

An  efficient  approach  to  preventing  transfer  of  acoustic  energy  from  one 
point  in  water  to  another  point  is  to  use  a  layer  of  material  which  is  more  compliant 
than  water.  Air  is  an  excellent  material  for  this  purpose,  as  shown  in  figure  3-1, 
where  the  insertion  loss  as  a  function  of  frequency  for  various  air  thicknesses  is 
plotted.  The  curve  shows  maximum  insertion  loss  at  \/4,  while  complete  transmis¬ 
sion  corresponds  to  multiples  of  a  half  wavelength  of  thickness. 

The  response  of  a  hydrophone  in  front  of  this  compliant  layer  must  also  be 
considered.  Figure  3-2  shows  the  response  of  a  hydrophone  1.5-inches  in  front  of  an 
air  layer.  A  severe  loss  in  response,  compared  to  free  field,  occurs  in  the  important 
frequency  band  of  1  to  3  kHz  due  to  the  180  degree  phase  shift  when  sound  is  reflected 
from  the  compliant  layer.  One  solution  to  this  problem  is  to  place  a  stiff  material, 
such  as  steel,  between  the  baffle  and  the  hydrophone.  When  steel  is  used  alone,  as 
shown  in  figure  3-3,  the  hydrophone  response  improves  considerably.  However,  the 
insertion  loss  of  steel  alone  is  low,  as  shown  in  figure  3-4.  The  combination  of  com¬ 
pliant  baffle  and  the  steel  has  only  a  small  affect  on  the  response  as  shown  in  figure 
3-5,  while  still  providing  a  high  insertion  loss. 

In  practice,  an  air  layer  by  itself  cannot  be  used  as  an  acoustic  baffle  for 
sonar  applications  because  of  the  pressure  requirements.  On  a  submarine  the  air 
layer  would  compress  at  a  moderate  depth  and  would  then  be  useless.  It  is  possible 
to  build  a  pressure-compensated  air  baffle  in  which  an  arragement  is  made  to  keep  the 
air  inside  at  the  same  pressure  as  the  water  outside;  however,  such  designs  are  usually 
complex  and  unreliable.  The  goal  is  to  build  a  baffle  which  is  statically  stiff  so  that  it 
can  withstand  the  hydrostatic  pressure  and  yet  dynamically  soft  so  that  it  provides  the 
required  insertion  loss. 

A  promising  approach  is  through  the  use  of  spring- mass  baffles  which  con¬ 
sist  of  a  sealed  assembly  of  plates  simulating  springs  with  the  outside  plates  behaving 
as  masses.  This  section  reviews  spring- mass  baffle  design  parameters  and  presents 
the  results  of  experimental  investigations  utilizing  full  scale  baffle  samples. 
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3.2  TRANSMISSION  LOSS  COMPUTATION 


The  configuration  of  an  uncompensated  spring- mass  baffle  developed  by 
Sperry  is  shown  in  figure  3-6.  The  assembly  utilizes  two  hardened  steel  inner  plates 
in  a  diaphragm  mode  as  springs  to  provide  high  compliance.  A  matrix  of  steel  pins, 
located  on  each  side  of  the  two  inner  plates  in  a  staggered  orientation,  transforms 
water  pressure  loads  from  the  outer  plates  into  concentrated  loads  on  the  inner  plates. 
This  arrangement  produces  diaphragm-type  inner  plate  deflections  in  the  form  of  con¬ 
tiguous  circles  in  each  inner  plate.  Each  circle  is  the  result  of  the  load  from  a  pin  on 
one  side  of  one  of  these  plates  being  resisted  by  the  pins  on  the  opposite  side  of  the  same 
plate.  The  pins,  held  in  place  by  square  rubber  pads,  fit  tightly  into  undersized  holes 
in  the  pads,  which  are  cemented  to  the  steel  plates.  The  top  outer  plate  is  supported 
from  the  bottom  outer  plate  by  a  perimeter  rubber  gasket  cemented  to  these  plates. 

Under  increasing  static  pressure,  the  inner  and  outer  plates  deflect;  the 
outer  plate  deflections  are  at  least  an  order  of  magnitude  smaller  than  the  inner  plate 
deflections,  both  because  the  outer  plates  receive  a  distributed  pressure  load  as  com¬ 
pared  to  the  concentrated  loads  from  the  pins  on  the  inner  plates,  and  because  the  outer 
plates  are  1/4-inch  thick  as  compared  to  the  1/8-inch  thick  inner  plates. 

The  lumped-parameter  model  of  the  spring- mass  baffle  is  shown  in  figure 
3-7(a).  In  this  figure,  the  transmitted  pressure  is  denoted  by  Pt  and  the  total  pres¬ 
sure  at  the  incident  face  is  denoted  by  Pq.  The  total  pressure  is  composed  of  the 
incident  pressure  Pj  and  the  reflected  puressure  PR;  as  follows: 


All  pressures  and  velocities  are  assumed  to  have  a  complex  exponential 
time  dependence.  Thus,  for  example,  the  actual  physical  incident-wave  pressure  is 
the  real  part  of  Pje-i^t,  where  u>  is  the  angular  frequency  of  the  incident  wave. 

Let  vj,  V2,  and  vg  denote  the  velocities  of  the  left-hand,  center  and  right- 
hand  masses  in  figure  3-7(a).  The  equations  of  motion  are 

-iwMVj  =  PQ  +  iK  (v2  -  Vj)/^ 

-iwmv2  =  iK  (vj  -  2v2  +  vg)/w 
-iwMv3  =  iK  (v2  -  v3)/w  -PT 

where  M  and  m  are  the  masses  per  unit  area  of  the  outer  and  inner  plates,  respectively, 
and  K  is  the  slope  of  the  baffle  pressure/displacement  curve.  By  interpreting  pressure 
as  voltage  and  velocity  as  current,  the  equations  of  motion  may  be  viewed  as  describing 
an  electrical  circuit  equivalent  in  behavior  to  the  mechanical  system. 
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1 

I 

The  equivalent  circuit  for  the  baffle  model  is  shown  in  figure  3-7(b).  The 
impedances  are  given  by 

Z..  =  -jtuM 
M  J 


Z  =  -iuim 
m  J 

Zr  =  jK/w 

The  effect  of  damping  in  the  springs  and  rubber  pads  of  the  baffle  may  be  recognized 
by  using  complex  values  of  K,  or  more  explicitly,  by 

ZK  =  L  +  ]K/oj 

where  L  is  the  effective  viscosity  of  the  spring/pad  combination. 

By  converting  the  interior  ^-section  of  figure  3-7(b)  into  an  equivalent 
T-section,  the  equivalent  circuit  for  the  baffle  is  converted  into  the  form  shown  in 
figure  3- 8(a),  where 


Z  .  =  Z„  +  Z  Zj{Z  +  2  Z„) 
A  M  m  K  m  K' 


JB 


<ZK>  /(Zm  +  2  ZK> 


When  the  baffle  is  immersed  in  water  and  subjected  to  an  incident  wave  of 
pressure  P],  the  pressure/velocity  relations  for  plane  progressive  waves  incident  upon 
and/or  radiated  from  the  baffle  surfaces  must  be  included.  At  the  right-hand  surface 
pure  radiation  exists,  for  which 

p  =  pt<>x/c 

where  c  is  the  speed  of  sound  in  water  and  x  is  measured  from  the  right-hand  surface. 


II 

i 

•  nip 


At  the  left-hand  surface,  both  incident  and  reflected  waves  exist,  so  that 


P  = 


iujx/c  „  -iu>x/c 
P  e  +  P  e 

I  R 


where  x  is  measured  from  the  left-hand  surface. 

A  water  particle  is  subject  to  the  general  law  of  motion 

dP/dx  =  -f)  dv/dt 

which  implies  that  P  =/>  cv  for  waves  moving  to  the  right  and  P  =  -p  cv  for  waves  moving 
to  the  left.  The  quantity  pc  is  the  characteristic  impedance  of  water.  Hence,  PT  may 
be  eliminated  from  the  equations  of  motion  by  the  substitution 


Pj  =  pcv  g  (20) 

l 
L 
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On  the  left-hand  surface,  the  following  exists 


PCV  =  P -  PDe-‘“x/c 

1  K 

so  that  by  setting  x  =  o  in  the  equations  for  the  left-hand  surface,  the  following  is 
obtained: 


PCV1  =  PI  -  PR 

from  which  is  obtained 

PQ  =  2Pj  -  pcvj  (21) 

Adjoining  (20)  and  (21)  to  the  equations  of  motion  shows  that  the  operation 
is  equivalent  to  the  circuit  shown  in  figure  3- 8(b),  where  Zc  is  the  characteristic 
impedance  of  water. 


by 


where 


Simple  circuit- analysis  techniques  show  that  the  transmission  loss  is  given 

TL  =  20  log  |  Pj/PT  | 


VPT  '  <ZA  +  ZC>  (ZA  +  2  ZB  +  ZC>/(2  ZB  ZC> 

Once  appropriate  values  of  M,  m,  K  and  L  have  been  chosen,  computer 
program  BAF,  which  is  listed  on  the  following  pages,  plots  TL  against  frequency. 

For  the  Sperry  baffle,  M  and  m  are  0.  00022  and  0.  000072  mass  units/in2, 
respectively.  K  varies  from  100  lb/in'5  at  shallow  depths  to  about  3000  lb/in3  at  700  ’ 
feet,  while  L  lies  in  the  range  0.2  to  0.  4  lb-sec/in3.  Figure  3-9  (Sheets  1  through  3) 
show  TL  versus  frequency  for  a  number  of  spring  constants  and  damping  values 
(A0M  =  M,  AIM  =  m,  AK  =  K,  AL  =  L). 


Experimental  data  on  a  spring- mass  type  of  baffle  has  shown  good  agree¬ 
ment  with  the  theoretical  calculations  over  the  low  frequency  band  (approximately  0.  5 
to  3  kHz).  At  higher  frequencies,  the  theory  predicts  larger  transmission  loss  than 
measured.  Also,  the  measurements  show  a  peak  and  a  null  in  the  transmission  loss 
curve  which  is  an  indication  that  the  baffle  thickness  becomes  equivalent  to  one-quarter 
and  one-half  wavelength  at  certain  frequencies.  When  this  occurs,  a  distributed  para¬ 
meter  model  is  required.  The  spring- mass  type  of  baffle  becomes  an  appreciable  part 
of  a  wavelength  at  even  relatively  low  frequencies  because  the  speed  of  sound  through 
such  a  baffle  is  much  lower  than  that  through  water.  Experimental  measurements 
indicate  the  velocity  through  the  baffle  is  only  one -sixth  that  of  water. 


LISTING  OF  COMPUTER  PROGRAM  TO  COMPUTE 
TRANSMISSION  LOSS  OF  SPRING-MASS  BAFFLES 


00100 
0  0110 
00  1P0 
00130 
0  0  1  40 
001  50 
001  60 
001  70 
00133 
00  1  90 
00200 
00210 
00220 
00230 
0  02  40 
00250 
00260 
0  02  7  3, 
00230 
00290 
00300 
0Q310 
00320 
00330 
00340 
00350+ 
00360 
00  3  70 
00380 
00390 
00  400* 
0041  0 
00420+ 
00430+ 
00  440  + 
00450+ 
00460 
00  473 
0 0  480 
00490 
00500 


PROGRAM  3AF( I NPUT, OUTPUT, TAP El  *  TAPE  6) 

DIMENSION  D3(70),FMTS(8,2),F0RM( 5),LINE(72),F(3,2> 

DATA  PI, RC/3. 1 41 59,5-84/ 

COMPLEX  ZI,Z2,Z3 
PRINT,  *WHAT  MODE,  MATE*, 

READ,  MODE 

DATA  F/ 1000., 200., 10000., 5., 5., 200./ 

REMIND  1  £  REWIND  6 

WRITE (6,500) 

530  FORMAT  (  4H - ) 

10READC1 ,510)  A0M, A  I M, NUMKL  £  I F  C  EOF, 1 )  110,15 
51QF0RMATC 4X2F10/4XI5) 

1 5D0  100  KL=I,NUMKL 
READ ( 1 ,520)  AK,AL 
520  FORMAT ( 4X2F 1 0 ) 

1=  1  $  FREQ=  F ( 1 ,M0DE) 

20T1  =  2  «  *P I +FREQ 

Z1  =  CMPLX<0.,-T1 *AIM) 

Z2=  CMPLX<AL,AK/T1 ) 

Z3  =  Z2/< Z1  +Z2  +  Z2) 

Z1 =  CMP LX  <  0 • , -  11 *A0M) +Z1  *Z3 
Z2=  Z2  *Z3 
G0T0C30, 40 )  MODE 

30  Z 1  =  Zl+RC  •  „ 

DBC I ) =  -  1000 . +A INTO  000 .5  +  20  . *AL0G1 0 ( CABS ( Zl* ( Z 1 +Z2+Z2 ) 
/(2 • *RC*Z2 ) ) ) ) 

GOTO  5Q 

40  DR  (  I  )  =  - 1 000.  +A I NT< 1000.5  +  20 . +AL0G1O CCA35C Z2/CZ1  +Z2) ) ) ) 
501=  1+1  5  FREQ  =  FREQ  +  FC2, MODE) 

IFCFREQ. LT.F(3, MODE)*. 5*F(2, MODE) )  GOTO  20 
WRITE  PLOT  HEADING 
DATA  FMTS/ 

50H( ////5X30HTRANSMI5S10N  LOSS  VS  FREQUENCY  >, 

30HC/4HCD3), I3,9I5,7H  (KHZ)  ), 

50 H( ////5X26K VELOCITY  GAIN  VS  FREQUENCY  ), 


3JH(/4H(03>, I3,8I5,6H  (HZ) 
ENC0DE( 50, 600, FORM)  (FMTSC I 
6OOF0RMATC  5A 1 0 ) 

WR ITEC  6, FORM) 

WR I TE ( 6, 610)  AOM, AIM,AK,AL 
6 1 0 FORMAT ( /7X, 4HAQM=,F8 .6,  7H 


)/ 

MODE), 1=1 


5) 


AIM=,F3.6, /7X4HAK  =i 


510  + 

F7.0,8H  AL 

= , F7 . 2 ) 

520 

ENCODE ( 30 , 60G , 

FORM)  ( FMTS ( I, MODE) 

,1=6,8) 

530 

1=  2 -MODE  S 

J=  1 90  +MODE - 1 80  S 

K  =  2  4*M0DE-23 

540 

WR I TE ( 6, FORM ) 

(L,L=l»J,K) 
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00550*FIND  D3  RANGE 

00560  DBMAX-  DBCI)  $  DBMIN=  DBCI) 

00570  1=  1  .5  +  CFC3*M0DE>-FC  WM0DE)>/FC2.M0DE> 

00580  00  60  J  =  2# I 

00590  D8MAX=  A MAX l C  DBMAX*  D3  C  J) ) 

00600  60 DBM I N=  A MINI C  D3MIN*D3C J) ) 

00  61 O*PL0T 

00  620.  Tl=  DBMAX 

00630  70LINEC 1 )=  6H 

00 640  IF CAM0DCT1 *5* ) •  NE  »Q • )  G0T0  80 
00650  ENC3DEC5*  620  >  L INE  >  T1 

00660  62OF0RMAT(F5.O> 

00670  8OD0  90  J=l#l 
00680  LINECJ+l )=  1H 

00690  IFCAM0DCT1 > 5* >*M0DC J+M0DE-2* 5> -EQ-O* )  LINECJ+l) 

00700  90IFCT1.E0.D3CJ>)  LINECJ+!>=  1H* 

00710  WR1 TEC 6* 640 )  LINE Cl )>CLINECL+1 ) . L  =  1 . I ) 

00720  640 FORMAT C A 6# 7 1 A 1 )  ■ 

00730  T1 =  T1 -1 . 

00740  IFCT1 .GE.D9MIN)  G0T0  70 

00750  1OOC0NTINUE 
00760  G0T0  10 

00770  1  10 MR  I TEC  6. 650) 

00780  650F0RMATC////4H - ) 

00790  ENDFILE  6  S  REWIND  6 

00800  END 
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3.3  EXPERIMENTAL  DATA 


Data  obtained  from  tests  of  a  spring- mass  baffle  and  associated  hydrophone 
array  are  presented  in  this  section.  The  array  tested  consisted  of  a  6  by  6  ft,  1-1  y2- 
inch  thick,  steel  plate  in  which  the  hydrophones  were  mounted,  a  1-1/2-inch  thick  rub¬ 
ber  sheet  bonded  to  the  front  of  the  plate,  and  the  spring-mass  baffles  placed  behind 
the  hydrophones.  The  test  array  contained  23  hydrophones  located  on  the  steel  plate 
as  shown  in  figure  3-10. 

The  following  tests  were  performed: 

•  Single  element  receiving  response 

•  Single  element  directivity  patterns 

•  Sum-difference  patterns 

•  Steered  beams. 

Test  were  conducted  at  depths  of  100,  400,  and  700  feet.  Of  the  large  number  of  test 
runs  made,  some  representative  samples  showing  the  baffle  performance  have  been 
selected.  Figures  3-11  through  3-16  show  the  single  element  directivity  patterns  of 
hydrophones  1  through  6  (see  figure  3-10),  which  are  located  across  the  array,  at 
2.  5  kHz.  It  is  interesting  to  observe  how  the  patterns  for  hydrophones  near  the  center 
of  the  array  (3  and  4)  are  symmetrical  while  those  at  the  edges  (1  and  6)  show  a  definite 
filling  in  of  the  pattern.  This  can  be  explained  by  noting  that  the  hydrophones  mounted 
near  the  array  edge  do  not  get  the  full  effect  of  the  baffle.  These  patterns  were  taken 
at  a  depth  of  400  feet  and  show  a  front-to-back  ratio  of  15  to  20  dB. 

Figures  3-17  through  3-19  show  broadside  beam  patterns  formed  by  sum¬ 
ming  the  outputs  of  hydrophones  1  through  6  at  1,  2.  5,  and  3.5  kHz.  The  theoretical 
beam  patterns  are  superimposed  on  the  measured  data.  Agreement  between  theory 
and  experiment  is  excellent  in  the  forward-looking  direction.  The  presence  of  the 
baffle  reduces  the  rear  response  by  a  factor  of  better  than  20  dB.  Depth  for  these 
tests  was  100  feet. 

Figures  3-20  through  3-24  show  steered  beam  patterns  at  30  degrees  and 
45  degrees  with  and  without  shading.  Hydrophones  1  through  6  were  steered  using  de¬ 
lay  lines;  shading  was  added  by  switching  in  resisitive  attenuators;  depth  was  700  feet. 
Theoretical  patterns,  calculated  assuming  free  field  hydrophones,  are  also  shown. 
Agreement  is  good  for  the  forward  looking  direction.  For  the  broadside  beams,  the 
baffle  attenuates  the  second  main  beam  by  more  than  20  dB. 

The  results  of  these  tests  indicate  that  the  spring-mass  baffle  is  capable 
of  effective  performance  over  a  wide  range  of  frequencies  and  hydrostatic  pressures. 
These  baffles  require  no  external  pressure  compensation  and  can  be  designed  to  be 
relatively  light  and  thin.  The  baffles  are  highly  compliant,  thus  providing  large 
transmission  loss  values. 
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Figure  3-1.  Insertion  Loss  Versus  Frequency  for  Various 
Air  Layer  Thicknesses 
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Figure  3-4,  Insertion  Loss  Versus  Frequency  of  Various 
Steel  Layer  Thicknesses 
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Figure  3-5.  Response  of  a  Hydrophone  1.  5  Inches  in  Front  of  a 

Combined  Air/Steel  Baffle 
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(a)  Lumped- Parameter  Model  of  Baffle 


(b)  Equivalent  Circuit  of  Baffle 


Figure  3-7.  Lumped  Parameter  Model  and  Equivalent 
Circuit  of  Spring-Mass  Baffle 


TRANSMISSION  LOSS  VS  FREQUENCY 

A0M=  .000220  AIM=  .000072 
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Figure  3-9.  Plots  of  Transmission  Loss 
Versus  Frequency  for  Various  Spring 
Constants  and  Damping  Values  (Sheet  2  of  3) 
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TRANSMISSION  LOSS  VS  FREQUENCY 
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Figure  3-9.  Plots  of  Transmission  Loss 
Versus  Frequency  for  Various  Spring 
Constants  and  Damping  Values  (Sheet  3  of  3) 
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Figure  3-16.  Directivity  Pattern  of  Hydrophone  6  at  2.  5  kHz 
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Figure  3-20.  30-Degree  Beam  Pattern,  Hydrophones  1  Through  6  at  2  5  kHz 

With  No  Shading 
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Figure  3-21.  30-Degree  Beam  Pattern,  Hydrophones  1  Through  6 
at  2.  5  kHz  with  20  dB  Shading 
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45-Degree  Beam  Pattern,  Hydrophones  1  Through  6  at  2.  5  kHz 
With  No  Shading 
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SECTION  4 


SONITE 


4.  I  INTRODUCTION 

,  T  .  So‘lite  ajj^hly  porous,  iiber-reinforced,  ceramic  material  manufactured 
by  the  Johns- Manville  Corporation.  It  was  originally  designed  as  a  thermal  insulation 
undei  the  trade  name  MIN-K,  According  to  the  manufacturer,  Sonite  has  been  speci¬ 
fically  formulated  as  a  low  acoustic  impedance  material  over  a  wide  pressure  ranee  for 
use  m  deep  submergence  applications.  ^ 

Dui  ing  the  manufacturing  process,  Sonite  is  compressed  to  yield  initial 
densities  ranging  from  20  to  35  pcf.  Under  static  pressure,  such  as  would  be  en¬ 
countered  in  an  ocean  environment,  both  the  density  of  and  speed  of  sound  in  Sonite 
increase.  In  order  to  give  repeatable  performance  under  pressure  changes  the 
material  must  be  prestressed  through  several  cycles  to  the  highest  pressure  it  is 
expected  to  encounter.  After  being  prestressed,  Sonite  exhibits  elastic  mechanical 
and  acoustical  properties  up  to  the  prestress  limit.  Its  acoustic  impedance  increases 
linearly  with  pressure  but  remains  low  compared  to  water.  Tests  by  Johns-Manville 
le  eronee  4)  have  shown  that  Sonite's  major  properties  are  unaffected  by  temperatures 
from  4  to  30  degrees  Centigrade.  1 

4.  2  THEORETICAL  CONSIDERATIONS 

The  major  component  of  Sonite  is  a  microscopic  particulate  material. 

Milan  amounts  of  fiber  and  phenolic  are  added  to  stabilize  its  structure  and  increase 
hand  liability.  Smce  thc  particulate  material  has  an  open  cell  structure,  the  cells  tend 
o  collapse  under  stress.  However,  if  the  Sonite  is  prestressed,  the  particles  will  not 
catastrophically  collapse  unless  the  prestress  is  exceeded.  Up  to  this  limit  the  indi¬ 
vidual  particles  undergo  reversible  deformation.  As  measured  in  reference  4  the 
properties  of  Sonite,  under  pressure  cycling,  become  stable  only  after  several  cycles. 
This  can  be  explained  by  realizing  that  during  each  stress  cycle,  the  Sonite  particles 
will  rearrange  themselves  until  an  equilibrium  position  is  reached.  The  reorientation 
process  is  believed  to  be  exponential  with  a  relaxation  constant  of  two  cycles 
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Expressions  for  Sonite  density  and  acoustic  velocity  as  a  function  of  pre¬ 
stress,  initial  density,  and  applied  stress  have  been  derived  (reference  5)  and  are 


where 


<>o 


1  -“ln  (VPi„>  -  B,„  r™  in  (P  /P.  )1" 

in  H  [_  nr  in  J 

Cin  +  ?  [°ln(VPi„>1m+  SI>1  1  '2 


p  =  density 

v  =  acoustic  propagation  velocity 
P0  =  initial  density 
Pm  =  prestress 

Pin  =  initial  stress  supported  by  the  sample 
P  =  applied  stress 
Bin  =  initial  compressive  stiffness 
B  =  compressive  stiffness 
Cin  =  initial  stiffness 

m,  n  =  experimental  parameters  within  the  system  related  to  detailed 
coupling  mechanism 

a,(i, y,S  =  proportionality  factors 


It  can  be  seen  that  the  properties  of  Sonite  can  be  changed  to  fit  a  particular  require¬ 
ment  by  applying  the  proper  prestress,  P  . 

The  effects  of  the  additional  components  used  in  the  manufacture  of  Sonite 
are  important.  At  low  pressures,  before  Sonite  has  firmly  compacted,  the  fiber  and 
phenolic  contribute  to  the  material's  mechanical  behavior.  At  moderate  pressures, 
because  the  load-bearing  particles  are  much  stronger,  the  effects  become  negligible. 
The  phenolic  also  acts  to  increase  the  shear  stiffness  of  Sonite  by  forming  weak  bonds 
between  the  microscopic  particles  locking  them  in  place.  Increasing  shear  stiffness, 
the  dominant  factor  in  Sonite,  forces  a  corresponding  increase  in  the  sound  velocity 
and  thus  raises  the  acoustic  impedance.  In  general,  anything  that  tends  to  bond  the 
particles  together  will  increase  velocity;  breaking  the  bonds  will  decrease  velocity. 

As  a  general  guide  to  the  use  of  Sonite,  in  order  to  avoid  problems 
associated  with  its  basic  structure,  it  is  recommended  that  the  maximum  applied 
stress  be  20-percent  less  than  the  prestress  and  at  least  five  cycles  of  prestress 
should  be  used  to  stabilize  the  material.  To  avoid  contamination,  which  can  form 
bonds  between  the  particles  and  degrade  Sonite 's  properties,  the  material  must  be 
carefully  handled  and  encapsulated. 
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4.3  EXPERIMENTAL  INVESTIGATION 


Semite  is  extremely  porous  and  must  be  protected  from  direct  contact  with 
water.  This  fact  complicates  the  problem  of  obtaining  test  data  which  reflects  the 
acoustic  properties  of  Sonite  alone. 

For  the  experiments  conducted  at  the  Sperry  Sonar  Test  Facility,  nine 
12  x  12  x  1-inch  panels  were  used.  These  panels  were  formed  into  one  three-foot 
square  structure  by  cementing  thin  (1/16  inch)  aluminum  panels  to  both  sides  and 
sealing  the  edges  with  cement.  At  the  depths  available  for  test  purposes  (12  feet)  no 
problems  with  leaks  were  encountered. 

The  test  configuration  is  shown  in  figure  4-1.  Two  series  of  runs  were 
made,  one  with  the  reference  hydrophone  1-5/8  inches  from  the  baffle  and  the  other 
at  a  distance  of  3-1  4  inches.  Frequency  response  data  from  1  to  20  kHz  were  taken 
with  the  baffle  and  hydrophone  located  as  shown  in  figure  4-1  and  with  the  baffle  and 
hydrophone  rotated  180  degrees  so  that  the  baffle  was  between  the  hydrophone  and 
projector.  Horizontal  patterns  at  various  frequencies  were  also  obtained  with  both 
hydrophone  spacings. 

Test  results  are  shown  in  figures  4-2  through  4-4.  Figures  4-2  and  4-3  are 
the  frequency  responses  for  hydrophone-to-baffle  spacings  of  1-5/8  inches  and  3-1/4 
inches,  respectively.  The  responses  at  180  degrees,  compared  to  free  field,  are  a 
measure  of  the  baffle  insertion  loss.  They  indicate  that  1  inch  of  Sonite  can  attenuate 
sound  in  the  1  to  20  kHz  region  by  at  least  10  dB.  However,  the  relatively  small  size 
of  die  test  baffle  causes  diffraction  to  limit  the  baffle  performance.  Figure  4-4  is  a 
comparison  between  measured  and  predicted  responses  at  a  hydrophone  spacing  of 
1-5  8  inches.  Agreement  is  good  over  the  1  to  10  kHz  frequency  range.  The  predicted 
results,  calculated  by  the  Sperry  multilayer  baffled- hydrophone  response  computer  pro¬ 
gram,  include  the  effect  of  the  two  1  16-inch  aluminum  panels  cemented  to  the  Sonite. 
However,  another  calculation,  with  a  baffle  composed  of  only  Sonite,  showed  that  this 
effect  is  negligible  due  to  the  small  thickness  of  the  aluminum.  The  acoustic  impedance 
of  Sonite  (20  pcf)  used  in  the  calculations  was  0. 132  that  of  water.  This  is  slightly 
lower  than  the  more  recent  value  of  0.  145  (after  10  cycles  of  0-1000  psi)  reported  by 
Johns- Manville  in  reference  4.  However,  the  Sonite  used  in  the  tests  had  not  been 
pressure  cycled  but  had  just  been  pressed  to  yield  a  density  of  20  pcf. 

Test  results  on  Sonite  to  date  are  limited  in  scope  but  indicate  its  useful¬ 
ness  as  a  baffle  material.  Its  light  weight,  coupled  with  an  apparent  ability  to  maintain 
a  low  acoustic  impedance  under  high  static  pressures,  makes  it  a  good  candidate  for 
future  investigation  and  experimentation. 
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,  ROTATION  (BAFFLE  a  HYDROPHONE) 

f  A  FOR  HORIZONTAL  PATTERNS 
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Figure  4-1.  Sonite  Test  Configuration 
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1 

I  Figure  4-2,  Relative  Hydrophone  Response  In  Presence 

of  Sonite  Baffle,  Spacing  1-5/8  Inches 
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Figure  4-3.  Relative  Hydrophone  Response  In  Presence 
of  Sonite  Baffle,  Spacing  3-1/4  Inches 
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Figure  4-4.  Comparative  Response  of  Sonite  Baffle 
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SECTION  5 


SUMMARY 


The  following  is  a  summary  of  the  results  of  the  investigations  made 
during  the  period  covered  in  this  report. 

•  A  baffle  composed  of  an  infinitely  long  cylinder,  with  walls  made  of 
layers  of  various  materials,  has  been  analyzed.  Results  of  the  analysis 
have  been  utilized  in  a  computer  program  which  calculates  the  response 
of  a  hydrophone  placed  anywhere  outside  the  baffle.  Limited  compari¬ 
sons  between  experimental  and  predicted  data  have  shown  extremely  good 
agreement,  although  more  work  remains  to  be  done. 

•  Studies  were  made  to  determine  if  a  uncompensated  spring-mass  type 
of  baffle  could  be  designed  to  meet  underwater  noise  reduction  require¬ 
ments.  Lumped  parameter  models  were  used  to  calculate  the  perform¬ 
ance  of  various  configurations.  Results  of  tests  performed  on  a  spring- 
mass  baffle  show  that  the  lumped  parameter  models  developed  in  this 
study  are  valid  for  predicting  the  performance  when  the  effective  thick¬ 
ness  is  small  compared  to  one  wavelength.  Excellent  agreement  between 
predicted  and  measured  data  was  obtained  at  the  low  frequencies.  The 
tests  further  showed  that  good  performance  can  be  obtained  with  these 
spring-mass  baffles  over  a  wide  range  of  frequencies  and  hydrostatic 
pressures. 

•  Sonite,  a  fiber-reinforced,  ceramic  material  manufactured  by  the  Johns- 
Manville  Corporation,  has  been  experimentally  studied  to  learn  about  its 
performance  as  a  low- impedance  baffle  element.  Tests  performed  at  the 
Sperry  Sonar  Test  Facility  on  a  Sonite  panel,  together  with  laboratory  data 
concerning  its  behavior  under  pressure,  indicate  that  this  material  shows 
promise  for  deep  submergence  noise  reduction  applications.  Further 
experiments,  especially  at  high  hydrostatic  pressures  and  using  large 
baffle  sizes,  are  need  to  fully  characterize  its  performance. 
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